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ABSTRACT 

We propose that the remarkable object Swift J1644+57, in which multiple 
recurring hard X-ray flares were seen over a span of several days, is a system in 
which a white dwarf was tidally disrupted by an intermediate mass black hole. 
Disruption of a white dwarf rather than a main sequence star offers a number 
of advantages in understanding the multiple, and short, timescales seen in the 
light curve of this system. In particular, the short internal dynamical timescale 
of a white dwarf offers a more natural way of understanding the short rise times 
(~ 100 s) observed. The relatively long intervals between flares (~ 5 x 10"^ s) may 
also be readily understood as the period between successive pericenter passages 
of the remnant white dwarf. In addition, the expected jet power is larger when 
a white dwarf is disrupted. If this model is correct, the black hole responsible 
must have mass <^ lO^M©. 

Subject headings: accretion, black holes, white dwarfs 
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Introduction 



On 28 March 2011, the Swif t Burst Alert Teles c ope (BAT) detected a most unusual 



object, Swift J164449.3+573451 flLevan et al 



2011 



Burrows et al. 



201l|). Although in 



many ways this object (whose name we abbreviate to Swift J1644+57) initially appeared 
to resemble a classical 7-ray burst, its light curve soon showed that it was quite different. 
Still bright more than 6 x 10^ s after the initial trigger (see Fig. [1]), in its initial activity 
it exhibited repeated extremely short timescale flares (see Fig. After holding roughly 
steady for ^ 700 s, the flare causing the BAT trigger rose a factor of 10 in flux over the 
next ^ 400 s. Less than 1000 s later, the flux had fallen by a factor of 20, and ~ 10^ s after 
that, although still detectable, the flux was only 0.5% of what it had been at the peak. 
Most surprisingly, there was a comparable flare ~ 50000 s later, similarly lasting for only 
~ 1000 s, and a third, slightly brighter than the first, ^ 60000 s after that. Like the first 
two, the third flare continued to show very short timescale variation. 

Several more brief flares of comparable brightness followed, likewise separated by 
periods of flux two orders of magnitude weaker. Gradually the duration of the flares 
stretched and their amplitudes diminished, until after ^ 2 x 10^ s the system began a long, 
gradual decline in flux in which it has remained bright enough to be detected for more than 
1.5 X 10'' s. Thus, there appear to be a number of characteristic timescales of variation, 
spread over quite a wide dynamic range: rise-times as short as ~ 100 s; flare durations 
~ 1000-10000 s; quiescent periods ~ 5 x 10'^ s long; and a total event duration of more than 
10^ s. After 10^ s the best power-law fit to the light curve is that the observed flux decays 
oc but its large variability and spectral changes mean that it might still be consistent 

with the expected t~^f^ bolometric decay. 

Still another surprise came when BAT data preceding the trigger were examined. 
Approximately 3 days earlier, there was a precursor whose peak flux was ~ 0.07 times the 
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Fig. 1.— Long-term Swift XRT light curves in the 3-10 keV band (blue: WT, red: PC). 
(Top) Linear in time representation of the first 300,000 s, illustrating the recurring brief 
flares that gradually widen. (Bottom) Logarithmic in time representation of the entire light 
curve as of 29 August 2011, five months after activity began. Both these light curves and 



those in the 



Repository (lEvans et al. 



bllowing figure we re generated using the graphing tools of the online Swift Data 



20071). 
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flux of tlie first flare (IBurrows et al. Il201ll ). Searches for earlier episodes of emission from 
the same source in archival data revealed only upper limits. 

Optical, near-infrared, and radio observations co mbine to show that the source lies 



within 150pc of the center of a galaxy aX z = 0.354 (ILevan et al. 



201 ih . The peak 



l uminosity associated with these flares (interpreted as isotropic) is then ~ 4 x lO'^* erg s ^ 



(IBurrows et al. 



1 X 10^® erg s ^ (IBurrows et al 



20111). Even at its fain test detected level, the (isotropic) luminosity is 



the flux by a factor ~ 100 is likely (IBurrows et al 



20 1 II) . Although relat ivistic beaming that could enhance 



20111 ). the beaming-corrected power is 



still very large: ~ 4 x 10 erg/s at the peak, and a total emitted energy ~ 5 x 10 erg 



wi thin the first 
by lBloom et al. 



xlOj s of the event if the XRT flux is ~ 1/3 of bolometric, as suggested 



(1201 ih . 



During flaring episodes, the observed uFi, is greatest between 10 and 100 keV; it 
appears to be strongly absorbed by interstellar gas below ~ 5 keV in the rest-frame. By the 
time the near-IR measurements were performed, 2-4 days after the event began, their fluxes 
were ~ 10~^ of the hard X-ra y flux at its peak, and ~ 10~^ of its flux during quiescent 



periods (IBurrows et al. 



20111 ). 



Within six weeks of its discovery, several models were proposed to explain it, most 
concentrating on a picture in which a main sequence star was :i (dallv disrupted by 
passing too close to a 10^-10'' Mr^ black hole ([ 



Cannizzo. Troja Sz Lodato 



2011 



Socrates 



31oom et al. 



2011 



Burrows et al. 112011 



201 ll ). Although these models differ from one 



another in some of their details, they share a common outline: A main sequence star is 
thoroughly disrupted as it passes near a black hole, and much of its mass is distributed into 
an accretion disk around the black hole. A powerful jet is then created, whose radiation, 
created by the synchro-Compton mechanism, dominates what we observe. In this picture, 
the rise time of the flares is interpreted as indicating (after allowance for relativistic effects) 
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Fig. 2. — Detailed time structure of the flares. (Top) Short timescale structure of the first 
flare. (Middle) Overview of the first flare. (Bottom) Short timescale structure in the flare 
~ 110,000 s after the BAT trigger. Note that the entire span of the data in the bottom 
panel is only 1500 s. 
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the size of individual knots in the jet, and these are related to the gravitational radius of 
the black hol e. The duration of th e entire event is thought to indicate either the expected 



( iRees 



1988 



t scaling (iBurrows et al 



Phinney 



201 if ) due to the orbital period distri bution of tidal streams 



19891 ) or the inflow time of the accretion disk (jSocratesll201ll ). 



Although plausible in many respects, this consensus model leaves a number of questions 
unanswered. The tidal disruption radius in units of the black hole gravitational radius is 



Rr/Rg = R.{k/ff'\MBH/M,f'/{GMBH/^) = 50{k/f)'/'M-'JlMT, 



T. 



where k is the apsidal motion constant (dete rmined by the star's radial density profile) and 



/ is its binding energy in units of GM^ / R^ (jPhinney 



19891 ). Here A^* is the mass of the 



star in solar units, and we have taken the approximate scaling that i?* oc M* on t 



ne mam 



seque nce. The ratio A;// is ~ 0.02 for radiative stars, but 0.3 for convective stars (jPhinney 



19891 ). The circular orbital period at the tidal radius is then proportional to A^* and 
independent of Mbh because the dynamical time at the tidal radius is always ~ (Gp*)"^/^, 



and on the main sequence the mean stellar density p* oc 



2. 



PoTh{Rj 



1.0 X lO^M^ s. 



(2) 



However, the actual orbital period of matter 



be considerably longer. As pointed out by 



captu red in tidal disruption is likely to 



ReesI (Il988[ ). if the disrupted star approached 



on a nearly-parabolic orbit, the tidal streams follow highly-eccentric elliptical orbits with a 
wide range of energies, and therefore of semi-major axes and orbital periods. Their mean 
energy is likely to be comparable to the star's self-gravitational binding energy because 
the gravitational force of the black hole does that much work expanding and disrupting 
the star. On the other hand, tidally-induced rotation and the gradient of the gravitational 
potential across the star can lead to some gas being trapped on orbits with semi-major 
axis as small as ~ Rt{Mbh / M^Y^^ , even while the typical stream's orbit is larger by a 
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factor {Mbh/M^Y'^. It follows that even the shortest orbital period is larger than that of a 
circular orbit at Rt by ~ [Mbh / M^Y^"^ , here a factor ~ 10"^, and the typical period likely 
another factor of ~ 10^ longer than that. Intersections between stream orbits could lead 
to conversion of orbital energy to heat, diminishing these orbital periods, but in no case 
would they become shorter than Porbl-Rr)- In fa ct, numerical simulation o f the disruption 



of a main sequence star by a 10 Mq black hole (lAyal. Livio fc Piran 



2000l ) shows a rather 



continuous accretion rate with a rise time of a few times 10^ s and an overall duration of a 
few times 10^ s, as expected from these analytic estimates. 

It is hard to reconcile these timescales to those seen in the lightcurve. The rise time in 
the consensus model is said to reflect the light-crossing time across the black hole's horizon, 
but it is not clear what dynamics link that quantity to triggering a flare, nor is there any 
natural explanation for the flare duration. The circular orbital period at Rt is comparable 
to the inter-flare interval, but the orbital period of the tidal streams following eccentric 
orbits, which is the timescale at which the t~^/^ decay begins, appears to be at least two 
orders of magnitude longer than the time at w hich the flares merged into a smo other 



lightcurve. These difficulties have led some (e.g.. ICannizzo. Troja fc Lodatd (1201 ll )) to pose 
special requirements on this model. They suggest that the pericenter distance must not be 
a great deal larger than the black hole's ISCO, so that inflow is largely dynamical. If so, 
the black hole mass is ~ IO^Mq. By the estimate of Equation [1], the pericenter distance 
must then be a small fraction of Rt- 

Another question is how exactly the accretion drives a jet. If, as is generally believed, 
jets associat ed with black holes are powered by some variant of the Blandfor d-Znaiek 



mechanism (IBlandford &: Znajek 



1977 



McKinney 



2005 



Hawley fc Krolik 



20061 ). substantial 



magnetic field must be attached to the black hole horizon. The Poynting luminosity in the 
jet is then ~ (pB^r'^c, where is a dimensionless quantity that depends on the field geometry 
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and increases with black hole spin parameter a/M, but is generally significantly less than 
unity. In this object, the field on the horizon must then be ~ 1 x l^^f^cj)"^/'^ L]/^ M^]j ^ G. 
When the system has a long lifespan (e.g., in AGN), the magneto-rotational instability stirs 
MHD turbulence in the accretion disk and builds the magnetic field; numerical simulations 
of MRI-driven MH D turbulence show that generically ~ 10 orbital periods are required 



to reach sat uration fIStone et al 



19961). Accre tion may also accumulate magnetic fiux on 



the horizon f Beckwith. Hawley fc Krolik 



20091 ): its build-up rate depends, of course, on 



the structure of the large-scale magnetic field. Here, the accretion fiow must plunge into 
the black hole on a dynamical time, and it is unclear whether such a strong field could be 
generated. 



Still another problem raised by this model is ho w to understand the 



Bloom et al. 



very large amplitude 



(120111 ). relativist ic jets 



and very rapidly-varying fiares. As recognized by 
in blazars behave very differently: Their characteristic fluctuation amplitudes are more 
typically factors of a few than factors of a few orders of magnitude, and their duty cycle at 
high flux is usually considerably greater than it is during the flrst several flaring episodes 
of this event. One could rephrase this question as, "What generates the extremely bright, 
compact, and short-lived knots that in this model are assumed to account for the flares?" 

In this paper, we propose an alternative version of the tidal disruption model that we 
believe holds some promise for supplying answers to all of these questions. We suggest 
that the star that is tidally disrupted is a white dwarf, not a main sequence star, and 
that it is not disrupted all at once, but instead loses pieces of itself in several passes 
before dissolving (see Fig. [3] for a schematic cartoon) an idea previously considere d in 



the context of white dwarfs on more n early circular orbits by 



Sesana et al. 



fl2008h and 



Zalamea. Menou fc Beloborodovl (120101 ). This suggestion is motivated by the fact that the 



fundamental timescale of a tidal disruption is dictated by the mean density of the star; 
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the greater density of a white dwarf makes it much easier to achieve the short timescales 
of this event. In both our picture and the main sequence star model, the accretion flow 
drives a jet, which produces the observed radiation; as we shall see, this, too, is favored 
quantitatively by the higher density of a white dwarf. The remainder of this paper will 
develop the consequences of these ideas. 



Note that tidal disruption of a white dwarf by a black hole has been also previously dis- 



cusse 


1 in the context of the possible nuclear ignition of the white dwarf ( 


Wilson & Mathews 


2004; 


Dearborn. Wilson & Mathews 


2005; 


Rosswoe. Ramirez-Ruiz & Hix 


2008 


). This 



requires a rather small black hole and a deep encounter, with pericenter a small fraction of 
Rt- Here we focus on larger black holes and more distant tidal disruptions that do not lead 
to such an explosion. 



2. Timescales 



Using the Using the mass-radius relations of iNauenberg] (119721 ). we find that a white 
dwarf is disrupted at 



(3) 



where Mj^ax — 1.4 is the maximal mass of a white dwarf in units of Mq. In this estimate. 



3) 



20091 ) and a binding energy 



we used an apsidal motion constant A; = 0.14 fISirotkin fc KimI 
factor / = 6/7, both appropriate to an = 3/2 polytrope. Note that we have changed our 
fiducial black hole mass from 10^ to IO^Mq because the higher densities of white dwarfs 
require smaller mass black holes in order to keep Rx/Rg > 1. In fact, white dwarf tidal 
disruption requires Mbh < 2 x IO^Mq, and somewhat less than that if the black hole 



rotates slowly or it is desired that Rt exceed the ISCO. The period of a circular orbit at 
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Fig. 3. — Top: A schematic description of our model, a white dwarf in a highly eccentric 
orbit that passes near a massive black hole. A significant fraction of the white dwarf mass 
is torn off around the tidal radius, Rt- The debris produces a small but massive accretion 
disk that supports a jet. Once the white dwarf moves away from the black hole the tidal 
disruption ceases, the disk is drained, and the jet dies out. Bottom: Different observed time 
scales and their relation to the schematic model. 
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the white dwarf tidal radius is (ignoring black hole spin) 



i"orb(i?T) ^ 67W;1 S, 



(4) 



three orders of magnitude shorter than for a disrupted main sequence star. In this regard, 
it is worth no ting that the white dwa rf mass distribution appears to be centered at 



O.6-O.7M0 flHansen fc Liebert 



20031). 



Matter torn from a white dwarf, like that taken from a main sequence star, will 
tr avel initially on a wide range of highly-eccentric orbits. Following the arguments 



of iReed (jl988l ). the most tightly-bound matter might be expected to have an orbital 

1/2 — S/"^ 

period ~ 900M^jj4A^* s. However, there are several potential mechanisms that could 
substantially shorten this period. If the black hole mass is not greatly smaller than its 
upper bound, the orbital pericenters will be i n the s trongly relativistic regime. There, 



as pointed out by 



Cannizzo. Lee fc GoodmanI fjl990l ). because relativistic effects cause 



differential pericenter precession, different stream orbits can intersect, so that dissipation 
transfers energy from orbital motion to heat and possibly photons. This process may be 
able to circularize the orbits on a timescale of a few orbital periods. When the black hole 
spins rapidly and its mass is close to the upper bound for tidal disruption, RT/Rg may be 
small enough that the orbital precession rate becomes comparable to the orbital frequency, 
enhancing this effect (we thank the anonymous referee for this suggestion). 

Multi-pass tidal break-ups may also lead directly to tighter orbits for tidally-disrupted 
matter. The orbital time of the white dwarf remnant is ~ 10 disk drainage times, so 
the only mass remaining in the disk by the time the remnant returns is the mass pushed 
outward absorbing the angular momentum of the accreted mass. If the mass remaining 
in the disk from the previous encounter is a fraction e of the newly-arriving mass, the 
momentum of the matter breaking off the star near Rt will be reduced by that fraction 
and its kinetic energy reduced by 2e. The semi-major axis of the resulting orbit for matter 
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with the mean incoming orbital energy (i.e., near zero) is then i?r/(4e), giving an orbital 
period Porb(-RT)(4e)~^/^. If e is as little as 0.05, the resulting period is only ~ lOPorbl-Rr)- 
Because there was a precursor to this event, it is possible that this sort of process may have 
occurred as early as the flare that caused the BAT trigger. 

As remarked before, internally-generated magnetic field requires ~ 10 orbital periods 
to reach saturation. When the period is as short as Porb(-Rr) for a white dwarf tidal 
disruption, the magnetic field can grow to full strength within ~ 100 s. Such rapid growth 
would then be consistent with explaining the rapid rise time of the flares. 

The ~ 1000-10000 s duration of individual flares might now be reinterpreted as the 
drainage time for the small, temporary accretion disk formed from captured material. Once 
the turbulence has reached saturation, this timescale tin ~ {PoTh/'^'^)(^~^{RT / if matter 
moves inward only by the action of internal disk stresses, where a ~ 0.1 is the ratio of 
vertically-integrated magnetic stress to vertically-integrated pressure, and the disk density 
scale height is H . Close to the ISCO, however, the mean stress overestimates the inflow 
time because stress fluctuations can remove enough angular momentum from fluid elements 
to send them all the way into the plunging regio n, permanently removing them from the 



disk. As a result. 



Krolik. Hawley fc Hirosd ( l2005l ) found that for R < SRjsco iii a disk with 



saturated MHD turbulence, tin ~ lOPorb(-R)- Allowing for some material on orbits with 
periods greater than Porbl-Rr), ^ 1000 s might then be a reasonable estimate for the total 
duration of a single flare, including both build-up of the magnetic field and inflow. 

If the initial pericenter of the white dwarf's orbit is ~ Rt, in its first passage it may 
lose only a fracti on of its mass (indeed the detection of a precursor flare ~ 3 days before 



the BAT trigger 



Burrows et al. 



(120111 )] suggests that the very first passage may have 
been slightly outside -Rt)- If that is the case, some of its mass is captured into orbiting 
streams by the black hole, but the rest of the star remains a coherent bound object, albeit 
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substantially stretched and distorted. It can then make several more passes through the 
region near R^, traversing a highly eccentric, but well-defined, orbit. 

If its initial approach was on a parabolic orbit, the white dwarf's net binding energy 
to the black hole is likely to be comparable to its initial self-gravitational binding energy 
because the energy to alter its structure is taken from the orbital energy. The period of its 
orbit after capture is then approximately 



-Porb ( Rc 



2n fM. 



'BH 



Rl 



3 \ 1/2 



6 X 10'^M;^Mbha s- 



(5) 



It is natural to identify this period with the interval between flares, as each time the 
star passes through the pericenter of its orbit it loses a fraction of its mass (see Fig. [3]). 
Moreover, because in each of these passes additional gravitational work is done stretching 
the remnant of the white dwarf, its orbit changes from pericenter passage to pericenter 
passage, explaining the irregularity of the inter-flare intervals. These structural changes 
may also alter the effective Rt from one passage to the next. 

When the white dwarf is flnally disrupted completely (presumably just before the 
flnal flare at ~ 1.7 x 10^ s after the initial BAT trigger), its remaining mass flnds itself 
spread over orbits with a wide range of binding energies and periods. The tidal streams 
with orbital periods longer than the remnant's orbital period return to the neighborhood 



of Rt only after a comparatively long time, extending the duration of the ev ent. 



distribution function of tidal-stream mass w ith orbital energy is roughly flat (iRees 
3 



the 



1988 



Phinne^ 



1989 



Lodato. King fc Pringle 



20091 ). the late-time accretion rate should decline 



oc t from the time at which the white dwarf is completely disrupted onward. 
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3. Accretion physics and driving the jet 

In the interior of a white dwarf, the electrons are highly degenerate. The 
equivalent temperature of the Fermi level Ep is 2 x lO^pg^'^ K, where we scale to a 
characteristic density of 10^ gm cm~^ because the mean density of a white dwarf is 



1.1 X 10^ M ^ M^l / (1 — 0.64:{M /MmaxY^^) gm cm ^ (again using the mass-radius re^ 



Nauenberg 



(jl972l )). On the other hand, typical interior temperatures are ~ 10^ K ( iHansen 



at ion o 



20041 ) ■ The initial phases of tidal disruption do not cha nge this degree of degene racy because 



adiabatic expansion leaves the ratio Ep/kT invariant (ILandau fc Lifshitz 



1980|) 



The flow does not remain adiabatic for long, however. For example, orbital precession 
leads to shocks. In these shocks, electron heating will initially be retarded by scattering- 
suppression due to their degeneracy, but there is no such constraint on the ions. The slowest 
ion-electron heating rate occurs when the ions are so much hotter than the electrons that 
the relative velocity is dominated by the ion thermal speed. Before allowing for electron 
degeneracy, the characteristic time for heat transfer by ion-electron Coulomb scattering in 
the disk formed by tidal disruption of a white dwarf is 

UonMat ~ 8 X 10-''{H/R)MBH,4T!ilM:%AM/M,)-' s. (6) 

Here we have set the Coulomb logarithm to 30 and assumed that, appropriate to a C/0 
composition, the mean ion mass is lArrip. The ratio AM/M^ is the fraction of the white 
dwarf mass deposited in a disk with radius Rp and scale height H . Degeneracy retards 
energy transfer to electrons because the crowded phase space partially suppresses collisions. 
Only electrons with initial momenta close enough to the Fermi level momentum that the 
momentum transfer per event Ap^e ~ ^e{vi) can lift them to the region of unoccupied 
states can participate. The fraction of electrons in the total population able to scatter is 
then 

/scatter ~ 3 K/2m,) ^fcT,/E^) _ Q^Q^^^TJ Epf'^ (7) 
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if Apie ^ Pf- The unadjusted heating time is so short that it is hard to imagine 
circumstances in which the degeneracy correction could make any difference on the 
timescales relevant to this situation. 

Thus, the electrons and ions can be expected to thermally equilibrate very rapidly, 
and at these densities and temperatures, the electrons (now no longer degenerate) will also 
rapidly thermally equilibrate with radiation. However, the radiative cooling time of the 
system is much longer than any of the relevant timescales: 



teooi ~ 1 X 10'\H/R)(AM/M,)M'J''M-^Jl s. 



(8) 



The radiation pressure is therefore effectively trapped within the material. Because 
shrinking the highly eccentric orbits to nearly circular requires dissipating an energy 
comparable to the orbital energy, the disk can therefore be expected to be geometrically 



thick, H/ R ^ 1. T 



accretion 



l is is t h e regime of photon- 



Begelman 



(119791 ): 



t rappin g associated with super-Eddington 



Abramowicz et al. 



( 1l988l ). Because the diffusion time is long 



compared to the inflow time, the radiation intensity distribution is far from steady-state, 
and the emergent luminosity is much less than the the rate at which heat is dissipated. 
Consequently, the ratio of thermal disk luminosity to rest-mass accretion rate is much less 
than the conventional ~ O.lMc^. 

The total pressure in the inner disk can be expected to be of order the electron density 
times the local virial temperature. 



dyne cm ^. 



(9) 



If the magneti c pressure on the event hori z on is limited by the inner disk pressure (the 



simulations of 



Beckwith. Hawley fc KrolikI ( l2009l ) suggest that it may be a factor of 3-4 



smaller), the expected field strength would be 

i?hor ^ 4 X lQ^\AM/M,f/'Ml"M-Jlll{H/Rr^ 



/2 



R 
10^ 



G. 



(10) 
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One might then predict a jet Poynting power as much as 

Ljet ~ 3 X 10^V(5L/87rpdisk)(AM/M,)7M,M^^^4 erg/s. (11) 

Even after allowing for a field rather less intense than the disk pressure, AM/M* < 1, and a 
small coefficient 0, it would seem that the jet power could easily reach the level seen (peak 
power after allowance for beaming of ~ 3 x 10^^ erg s~^). Note, also, that Lj^t oc M^jj, so 
models requiring larger black holes tend to generate weaker jets. As pointed out by our 
anonymous referee, the smalle r (in Rg terms') disks characteristic of white dwarf disruptions 



also minimize Compton drag (jPhinney 



19871) 



Numerical general relativistic MHD simulations of accretion have shown that the 
luminosity of the jet can also be estimated in terms of a n effective "efficiency" per unit rest - 



mass accreted that is a function of the black hole spin (IMcKinney 



2005 



Hawley fc Kroli 



20061 ). Using the expression found in the latter reference, one might predict 

Ljet ~ 4 X 10^9 [1 - \a/M\r^ {AM/M,)Ml erg/s (12) 

if the inflow time is ^ lOPorb(-Rr)- For a black hole with spin parameter a/M ^ 0.9, this 
estimate agrees with the previous one evaluated for the flducial parameters. That it should 
do so is no coincidence — as assumed in the previous estimate, the magnetic fleld intensity 
on the horizon is comparable to the inner disk pressure in these simulations. However, it 
should be borne in mind that the energy for this jet is actually drawn from the reducible 
mass (the rotational kinetic energy) of the black hole, not the accretion flow. The function 
of the accretion is solely to sustain a strong magnetic fleld on the black hole horizon. 

Although nominally independent of black hole mass, in fact this second estimate has 
an implied dependence through the bound on black hole mass placed by the density of 
the disrupted star. In rough terms, the second jet luminosity estimate simply mirrors the 
accretion rate, which is oc AMQ{Rt). In both the main sequence star and white dwarf 
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models, the amount of mass placed in orbit is ~ Mq. Where they differ is exactly the point 
we have emphasized in the context of the lightcurve's timescales: VL{Rt) ~ (Gp*)^/^, which 
is ~ 10'^ times larger for white dwarfs than for main sequence stars. 

When the disk is drained, the inner disk pressure falls. If that permits the flux on 
the horizon to expand, so that the magnetic field there becomes weaker, the jet would be 
correspondingly diminished. Without continuing accretion — that is, between episodes of 
tidal capture — radiation from the jet would be much reduced. This is, of course, what one 
would expect in those intervals when the remnant of the white dwarf is out near apocenter. 

4. Summary 

We propose that the remarkable event known as Swift J1644+57 is more likely the 
tidal disruption of a white dwarf than a main sequence star. The fact that white dwarfs are 
typically ~ 10^ times denser than main sequence stars makes it much easier to understand 
the very short timescales characteristic of this object's lightcurve. Postulating that the 
initial encounter was not close enough to disrupt the white dwarf completely explains the 
remarkable flares seen over the event's first few days (see Fig. [3]). The ~ 100 s fluctuations 
during the flares may then be identified with the inflow timescale from radii ~ Rt-, which in 
this case is likely only ~ IQRg. This timescale is imprinted on the lightcurve if the radiation 
emerges from the jet when it has traveled less than a few hundred seconds (in the observer's 
frame) from the black hole, a distance equivalent to ~ lO'^-lO'^-R^. The disk drainage time 
is longer because the disk can be expected to spread to radii somewhat larger than Rt'-, this 
accounts for the flare durations, ~ 10^-10^ s. The interflare time ~ 5 x 10^ is the orbital 
period of the white dwarf remnant. Because the remnant orbital period is rather longer 
than the drainage time, there is an extended period after most of the disk mass has been 
accreted, but before it is refilled by the next pericenter passage of the white dwarf, when 



- 19 - 



the accretion rate is low, and the jet is therefore weak. 



The higher densities of white dwarfs also explain both why the observed radiation is 
dominated by a jet and the jet's high luminosity. The very high surface density of the 
disk thoroughly traps thermal photons, while the associated high pressure can support a 
stronger magnetic field on the black hole. Although the jet luminosity associated with a 
main sequence star event is nominally just enough to supply the observed radiated power 
(after beaming corrections), that may still be inadequate. The numbers quoted refer only to 
luminosity in observed bands; it is possible there is additional luminosity elsewhere in the 
spectrum (e. g., between 100 keV and 1 00 MeV). More importantly, in most jet radiation 



models (e.g., ICelotti &: Ghisellinil ( l2008l )). the photon luminosity is only ~ 1% - 10% of the 
kinetic power; the jet powers estimated on the basis of the disk pressure or accretion rate 
refer to the total, so that Ljct must be significantly greater than the observed luminosity. 

Although the space density of white dwarfs (^ 3 x 10"'^ pc~^ in the Solar neighborhood: 
Rowell fc Hamblyi (120111 )) is very similar to the space d ensity of solar-mass main sequence 



stars (~ 3.5 x 10 ^ pc' 



-3. 



Reid. Gizis fc Hawleyl ( 120021 )). their effective cross section for 



coming close enough to a black hole to be tidally dis rupted is smalle r by t 



le ratio of their 



-2/3 



Burrows et al. 



(120 111 ) estimated the 



tidal disruption radii, ~ 2 x 10"^(A^^,^„d-^*,: 
total number of solar-mass main sequence star disruptions within the volume detectable by 
Swift to be ~ 10^ yr~^, suggesting that the rate of white dwarf tidal disruptions should be 
~ 10 yr~^. If they are all relativistically beamed so that we see only ~ 1% of all events, 
the observable rate falls to ~ 0.1 yr^^, and is re duced further by a fac tor ~ 10 to account 



for BAT's field of view and observing efficiency (IBurrows et al. 



201ll ). Thus, we might 



expect BAT to detect perhaps ~ 0.01 yr ^, rather than the actual ~ 1/6 yr ^. On the 



other hand, the rates estimated by 



Burrows et al. 



(120111 ) (and similarly by 



Zauderer et al. 



(120 111 )) were based on theoretical predictions that are an order of magnitude smaller than 
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some empirically-based rate estimates fiGezari et al. Il2008t iMaksym. Ulmer fc Eracleous 



20101 ). Given the many uncertainties (and small-number statistics), we believe the rate of 
white dwarf disruptions is consistent with the Swift detection rate. Moreover, a large part 
of the contrast between our predicted rate and the rate expected for main sequence stars 
would be removed if the r equired pericenter distance for a m ain sequence tidal disruption is 



0.1-Rt, as suggested by 



Cannizzo. Troja fc Lodatd (120111 ) 



We close with a final significant contrast between white dwarf and main sequence star 
models for this event. The black holes most effective at tidal disruption of white dwarfs are 
smaller by the ratio (p*,ms/p*, ^wd)^ ~ 10 ^. Although there has long been good evidence 
for black holes of ~ IO^Mq or more in the central regions of Galaxies (beginning with the 
Milky Way) and some what shakie r evide nce for black holes in galactic nuclei with masses 



^ IO^Mq (collected in 



Xiao et al. 



(120111 ). these masses are all based on the assumption that 
gravitational dynamics dominate broad-line gas motions, and in most cases use only scaling 
arguments to estimate the broad-line region's size), this event represents the first indication 
of a black h ole in the ^ lOMO^M ^ mass range. On the basis of the M^jif-bulge luminosity 



correlation, 



Burrows et al. 



(120111 ) estimated that this galaxy should have a nuclear black 
hole with mass ^ 2 x 1O^M0. If we are correct in our suggestion that this event was due 
to disruption of a white dwarf rather than a main sequence star, the correlation estimate 
may be too large by a factor of 100 or more in this galaxy. Alternatively, the black hole 
responsible for this disruption could be a second, smaller black hole orbiting a larger one. 



whose mass is closer to the correlation estimate. 
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